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ABSTRACT 

The most extended, closed magnetic loops inferred on T Tauri stars confine hot, X- 
ray emitting plasma at distances from the stellar surface beyond the the X-ray bright 
corona and closed large-scale field, distances comparable to the corotation radius. 
Mechanical equilibrium models have shown that dense condensations, or "slingshot 
prominences" , can rise to great heights due to their density and temperatures cooler 
than their environs. On T Tauri stars, however, we detect plasma at temperatures 
hotter than the ambient coronal temperature. By previous model results, these loops 
should not reach the inferred heights of tens of stellar radii where they likely no 
longer have the support of the external field against magnetic tension. In this work, 
we consider the effects of a stellar wind and show that indeed, hot loops that are 
negatively buoyant can attain a mechanical equilibrium at heights above the typical 
extent of the closed corona and the corotation radius. 
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1 INTRODUCTION 

Observations of T Tauri Stars (TTS) have revealed the pres- 
ence of magnetic structures much larger than those seen on 
our Sun. Indeed, analysis of these structures has invoked 
solar flare -based mode l s: to infer the length scales of confin- 
ing loops. [Reale et all (| 19971 ) developed the uniform cooling 
loop (UCL) model, linking solar flare X-ray decay slopes to 
spatially resolved lengths of the loops confining the emitting 
gas. The UCL model was applied to data from the Ch andra 
Orion Ultradeep Project (COUP. iGetman et al]|2005l ). and 
it was discovered that evidently the post-flare loops of TTS 
can r each up to multiple stellar radii in extent (|Favata et alj 
120051 . results summarised in Figure [T]). This result is surpris- 
ing: T Tauri coronae are typically believed to be compact, as 
evidenced by the common de tection of rotationally modu- 
lated coronal X-ray emission (Fla ccomio et al]|2005h . X-ray 
bright regions are generally thought to be confined within 
the corona, with loops generally not reaching more than a 
stellar radius in height. 

In the mechanica l equ ilibrium model of 
IJardine fc van Ballegooiier] ( 2 005), prominences form 
after reconnection occurs in the open field, and it is a 
straightforward matter to show that cool, dense loops can 
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reach heights well beyond corotation (the radius at which 
the orbital velocity equals the stellar rotational velocity) 
and the closed stellar coronal field. Hot loops such as 
those inferred by the COUP, however, appear to need the 
support of the stellar field, and thus should be limited to 
heights within the closed coronal magnetic field: a fraction 
of the heights inferred by UCL models. It was initially 
suggested that perhaps these hot, large-scale loops are 
anchored to circumstellar disc material for stability. Testing 
this idea, followup studies of the COUP objects that 
searched for discs were unable to find conclusive evidence 
for dusty disc material within reach of the loops (e.g., 
IGetman etafl l200Sl : [Aarnio et al]|20ld h Furthermore, the 
apparent anti-correlation of large-scale loops and close-in 
disc material suggests it is the absence of inner disc material 
that allows the loops to reach greater extents than they 
would were close-in disc material present (as suggested by 
IJardine et al. 2006, close-in circumstellar discs could strip 
away these stars' outer coronae). Potentially, this indicates 
some missing physics in earlier models, and it raises the 
question of how such large loops remain stable for the 
duration of the X-ray flare decay phase (in some cases, 
multiple rotation periods). 

The potential ramifications of highly energetic events 
on stellar and circumstellar evolution are far-reaching: large- 
scale magnetic structures, in their formation and eventual 
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Figure 1. A summary of the results of iFavata et al. | J2005T ): in- 
ferred loop height as a function of measured loop temperature. 
Vertical lines connect the loop height to that star's corotation 
radius. We have scaled the figure for ease of viewing; the tallest 
loop (its vertical corotation line shown) is ~63 R* in height. Filled 
points denote cases in which any parameters in calculation of the 
corotation radius were assumed (i.e. , any values of R* , M, , or P ro t 
not drawn from surveys of the ONC; fiducial values were used as 
discussed in section 13.111 . The grey horizontal line at 3 R* indi- 
cates the location of the source surface (the extent of the closed 
coronal magnetic field) in our models. There is a pileup of points 
at 270 MK because, to be conservative, the authors adopted a 
maximum best-fit temperature and applied it to all cases above a 
given threshold (beyond this temperature, the fits to their X-ray 
data are statistically equivalent, and a higher temperature simply 
yields a shorter loop). 

disruption, could regulate stellar angular momentum evolu- 
tion as well as circumstellar disc evolution, and thus planet 
formation. The combination of the high frequency with 
which large-scale magnetic structures have been inferred and 
the great heights at which they form above the stellar surface 
indicates circumstances are ideal for the generation of sub- 
stantial torques opposing stellar rotation. Large-scale mag- 
netic reconnection, in addition to shedding mass and angular 
momentum, also bombards circumstellar material with high 
energy photons and particles. This could potentially play a 
pivotal role in planet formation by stim ulating grain growth 
(cf. fl ash-heated chondrule formation, iMiura fc Nakamotd 
l2007h . High energy radiation and particle fluxes can infLu- 
ence disc chemistry (cf., solar nebu la isotopic abundances , 
iPeigelson et~aT1l2002l ) and structure: Idassgold et all (|l997l ) 
showed that hard stellar X-rays could potentially ionise the 
disc sufficiently for accretion to occur via magnetohydrody- 
namic instability. The large heights inferred for these loops 
would put X-ray emitting plasma closer to the disc, reduc- 
ing the column depth through which the hard X-rays would 
have to travel before ionising disc material. 

Given these potential impacts of large scale magnetic 
loops (and their destabilization) on early star and disc evolu- 



tion, we aim to ascertain whether the existence of hot, large, 
magnetic loops, as inferred from the analysis of the brightest, 
most energetic events in X-ray surveys of TTS, is physically 
plausible. To do this, we address an evident buoyancy de- 
fecit in the previous models by including a stellar wind in 
the model physics. We do not distinguish between the single 
loop UCL model or the possibility of multiple loop events; 
we only seek to demonstrate that it is possible to find hot 
loops beyond the typical "quiescent" X-ray emitting corona, 
and in some cases, beyond corotation. 



2 MECHANICAL EQUILIBRIUM MODEL 

Cool, extended prominences have been observed in abun- 
dance on the young rapid rotator AB Dor, rang- 
ing in height from 2-5 R*, well above the coro- 
tatio n radius of 1.7 R* (Col lier Cameron fc Robinson! 
1989), and slingshot prominences much like those 
seen on the Sun have been obser ved on TTS (e.g., 
iMassi et all 120081 ; ISkellv et all l2008h. From a theoreti- 
cal p e rspective, other groups | Jardine fc Collier Cameror] 



pe 

Il99ll ; IJardine fc van Ballegooiienl 120051 . hereafter JCC91 
and JvB05, respectively) have assessed the stability of cool 
loops, determining where in parameter space mechanical 
equilibria can be found. Loops found beyond the closed coro- 
nal field and possibly even beyond corotation are likely em- 
bedded in the stellar wind (these loops could form via re- 
connection of open field lines in the wind; see JvB05, Figure 
2). 

This work builds upon the prescription of JCC91 and 
JvB05, save our addition of a stellar wind (Section 12. 2p . We 
list below the four major forces our model accounts for (and 
the parameters upon which each depends) that dictate the 
height and shape of a magnetic loop: 

(i) Gravity (M», R», w), 

(ii) Buoyancy (p, g, T, dp/dy), 

(iii) Magnetic pressure gradient (assumed field structure), 

(iv) Magnetic tension (footpoint separation, external field 
strength and structure). 

A magnetic loop is in equilibrium when the pressure gradi- 
ents internal and external to the loop, magnetic tension, and 
gravity are evenly balanced: the loop is stable, neither ex- 
panding nor contracting. We describe the general mechanical 
equilibrium model, and the addition of a stellar wind term. 
In Section |31 the case of a simple TTS is demonstrated. 



2.1 Equilibrium Conditions 

We incorporate a number of fairly standard, simplify- 
ing assumptions in our 2D model. First, the arcade in 
which the loop is embed ded is not twisted (i.e. no shear; 
iBrowning fc Priestl [l983 . found, in the solar case, that in- 
cluding twist in the flux tube only serves to reduce the max- 
imum loop footpoint separation) . Secondly, the external field 
is potential (i.e. current free) and closed up to a height y s ; 
beyond this point, the source surface, the external field is 
open. The source surface sets, then, the maximum height 
quiescent coronal loops can reach. The magnetic loop itself 
is assumed to be narrow, its width less than the length scale 
of the external field, and it is isothermal. Finally, the loop is 
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Figure 2. A graphical representation of the form of the external 
field (see equation [2j . Our chosen source surface height, where 
the magnitude of the field is zero, is 3.0 R»; beyond this height, 
the field is open. Vertical lines indicate the width of the arcade 
(\x\ = ir/k, where here we choose k = 3 so the loop covers 60° in 
longitude about the stellar equator). 



treated as non-disruptive of its surroundings. The equation 
of motion for plasma in the corona is: 







PS, 



(1) 



where the first term expresses the plasma and magnetic pres- 
sure gradients, the second term represents the magnetic ten- 
sion, and the final term is the balance of gravitational and 
centrifugal forces. In all the following derivation, we refer to 
quantities inside the loop as "internal" , and coronal proper- 
ties as "external" (variable subscripts i and e, respectively). 

We choose a simple external field which is potential and 
two dimensional: 



Bex + lB ey = Bo (e^ 2 + I 



-2fcy a sl/2 



(2) 



where z = x+\y (see Figure [2] illustration of field) and Bo 
is the field strength at the stellar surface. Rewriting the 
Lorentz force in terms of the current density, j, we have, 



+ _( B .V)B=jxB. 



(3) 



Since the Lorentz force is zero in the direction of the field 
(hence j x B =0), we can combine equations ([T} and ((3} to 
obtain an expression for pressure balance along the loop: 

Ts =P9s > (4) 

where g s = g ■ B/|B| is the component of gravity along the 
loop. In a Cartesian coordinate system where g = (0,g), this 
reduces to: 



dp 

Ty = pg(y). 



(5) 



We assume the plasma along a given loop is in hydro- 
static equilibrium, so we solve for it as: 



p = poexp 



^ k B T 



gdy 



(6) 



where po is the loop base pressure (i.e., where y=0). 

We note that pressure balance across the loop requires: 



Bt = B; + 8n(p e - Pi ). 



(l) 



Following the derivation and boundary conditions of 
JvB05, we consider gradients normal to the field by tak- 
ing the scalar product of n and both sides of equation @, 
arriving at an expression that governs the shape of the loop, 



d \ Bf 



y 8x dy J 2 



1 + {y'Y 



B 



(8) 



The loop shape (given by y(x)) is therefore determined 
by the requirement that the Lorentz force is zero and hence 
that the magnetic pressure gradients are balanced by the 
magnetic tension. At the loop summit, magnetic tension 
acts directly downwards, so we must have dB 2 /dy < 0. The 
behaviour of Bf, however, is determined through pressure 
balance (equation by the nature of the external magnetic 
field and the variation of the internal and external plasma 
pressures. The requirement that 8B 2 /dy < reduces to: 



1 dB 2 e 
8tt dy 



< (Pi ~ Pe)i 



(9) 



From equation we begin to see some necessary phys- 
ical conditions for stable loops. The maximum loop height is 
constrained by how buoyancy (i.e., the gradient of the pres- 
sure difference, d(p e — Pi)/dy) depends on height. In the 
absence of a wind, the gas pressure external to the loop falls 
with increasing height until a turning point is reached at 
corotation: at this radius, centrifugal forces dominate over 
gravitational, and the buoyancy changes sign. If equation 
((9| is satisfied, the loop can be supported by the external 
field; if not satisfied, the loop will lose equilibrium. 

Meeting these conditions becomes a problem in the T 
Tauri case, as densities measured in hot loops are greater 
than ambient coronal densities, the internal temperatures 
are likely greater than the external temperatures, yet the in- 
ferred heights put them beyond corotation. By our previous 
understanding of the physics, this combination of parame- 
ters would not result in a mechanically stable loop within a 
hydrostatic corona. 

For cool prominences, equilibria existed beyond coro- 
tation only for negatively buoyant loops; in addition to the 
loop density being greater than its surroundings, the plasma 
/3 (ratio of gas to magnetic pressure) is sufficiently large, but 
less than unity. We can adjust the plasma /3 in the previous 
models and find solutions for tall, hot loops, but the value of 
P required for solutions beyond corotation is unphysical; for 
P values consistent with measured parameters on TTS, hot 
loops are positively buoyant to large distances. To provide 
additional upward force to support the hot loop, we add 
a stellar wind external to the loop. In this case, hot loops 
can reach a greater range of heights because the wind term 
extends the range over which buoyancy is negative. 
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2.2 Including a Stellar Wind 

In order to maintain stability of very extended, hot loops, 
the external plasma pressure must decrease approaching and 
passing corotation in order to allow the loop to maintain neg- 
ative buoyancy at large heights. Far enough from the star, 
the external magnetic pressure becomes negligible, and so 
the internal field strength (and thus the loop shape) is de- 
termined by the plasma pressure difference. Now, we include 
an additional term to include a lowering of the pressure due 
to the flow of a stellar wind along the open external field 
in which the magnetic loop is embedded; this will provide 
more upwards force on the loop, potentially allowing greater 
heights to be obtained by loops with a wider range of param- 
eters (consistent with inferred parameters of TTS magnetic 
loops) . 

We have added t o the basic formulation of JvB05 a 
simple, Parker wind (|Parkedll964h . Assuming an inviscid 
parcel of gas is corotating with a star, the total force on 
that parcel is: 



-Vp + pg+j x B. 



(10) 



As this is an expression of Newton's second law, it can be 
rewritten as follows: 

F = P^, (U) 
where is the substantive derivative, 

TT^l + ^.sj). (12) 

We will assume here both a potential field and a steady flow, 
so the momentum equation along the direction of the loop 
(i.e., parallel to the direction of the field: § = B/|B|) is: 



ps ■ (v ■ V)v = — s ■ Vp + pg ■ 



(13) 



Under the assumption of an isothermal gas, density 
and pressure are related by the local sound speed of the 
plasma, p — 2%. Conservation of mass requires that the 
flow across a given cross sectional area remain constant 
(i.e., pvA — constant,) and similarly the conservation of 
magnetic flux dictates no sinks or sources of magnetic field 
(BA — constant). 

Bearing these relationships in mind, the momentum 
equation along the loop can be re-written as: 

,,2 ^ 



d_ 

ds 



-<£ dp , . 
— + g-s. 

p OS 



(14) 



Integrating equation (|14[) from the loop base upward, 
we obtain the following expression relating the external field 
and pressure, and M = ^ the initial sonic Mach number of 
the flow: 

= \n(p e ) + Ul 2 - hi 2 - 1 / gdy. (15) 



The external pressure and field, p e and B e , are scaled by 
their base values (po and Bo, respectively). 

If the initial Mach number is 0, the momentum equation 
reduces to the hydrostatic case (equation [6}. The only phys- 
ical solution of equation (|15|l is the transonic solution, which 
gives a slow, sub-sonic speed close to the star, passes through 
a sonic point (where the flow velocity equals the local sound 
speed), and is supersonic at greater distances from the star 



(see Figure 1 of |Parkedll965l ). To illustrate how including a 
stellar wind can impact the maximum possible loop height, 
we consider the specific case of a TTS below. 



3 SPECIFIC CASE: A T TAURI STAR 

We have applied the mechanical equilibrium model plus 
wind to a generic, simple T Tauri star. The basic outline 
of the model and analyses over robust par ameter spaces are 
descr ibed fully in previous work (JCC91; lUnruh fc Jardina 
Il997l. JvB05); here we give full attention to the new addi- 
tion of a stellar wind term and its effects on the physics of 
the model. First, we discuss possible equilibria for hot loops 
embedded within the closed corona; second, we assess maxi- 
mum heights for loops extending beyond the source surface. 

3.1 Fiducial T Tauri Parameters 

As input parameters to the model, we require stellar mass, 
radius, and rotation period; we have adopted representa- 
tive mea n values from the Orion Nebula Cluster from the 
study of iHillenbrandl (|l997h . The mean mass of members 
(with 50% membership probabilities or greater) is 0.55 Mq, 
while the mean radius is 1.7 R,0. Our fiducial TTS rotation 
period of 5.3 days is derived from the assembled rotation 
period data of iGetman et al.l (|2005l ). We adopt a wind tem- 
perature of 2 MK, co nsistent with T Tauri wind observations 
l|Dupree et al.| [2005). Internal loop temperatures have been 
taken from the range of values measured f or the COUP "su- 
perflares" inferred by iFavata et alj (|2005h : these vary from 
~50-200 MK (see Figure [TJ. We set the initial Mach number 
to 0.05. For the extent of the closed coronal field, we show 
our results for y s = 3 R* , which is within the corotation ra- 
dius (5.2 R»). It is unclear for TTS how far from the stellar 
surface the closed corona may extend, but as we note in Sec- 
tion [4] our results do not depend strongly on the choice of 
source surface height (see Section [3.31 for further discussion 

ofy.)- 



3.2 Pressure and Buoyancy 

For the hydrostatic case, the functional form of the exter- 
nal plasma pressure is exponentially decreasing close to the 
star, with the argument of the exponential function depen- 
dent upon gravitational and centrifugal forces. In the case 
of a T Tauri star specifically, given the same mass and 
rotation period as a main-sequence counterpart but twice 
the radius, the effective gravity is approximately four times 
weaker, thereby extending the distances to physically impor- 
tant heights. The corotation radius for the TTS is a factor 
of two closer, and the difference in pressure (i.e., p e — Pi) 
remains below the base pressure difference to a distance ap- 
proximately three times closer than the main sequence star. 
The corotation radius is important because it is where pres- 
sure difference and the buoyancy change sign. 

We choose a base pressure and temperature ratio, and 
with these we solve equation (|15[) for the equilibrium shape 
of the loop. In Figure[3]we show, for the parameters outlined 
above, the internal pressure, pressure difference (p e — Pi), 
and buoyancy (d(p e — Pi)/dy) as a function of height for 
various temperature ratios T e /Tt. It is important to note 
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Figure 3. Pressure, pressure difference, and buoyancy as a func- 
tion of height (upper, middle, and lower panels, respectively). The 
external plasma temperature is 2 MK, and temperature ratios 
(T e /Ti) shown arc 1.5 (solid line), 1.0 (dotted line), 0.75 (dashed 
line), and 0.5 (dot-dashed line). The corotation radius is located 
at 5.2 R* (grey vertical lines); at 20.2 R, the internal pressure 
equals its base value, and the pressure difference changes sign. 



the heights at which the loop's buoyancy is negative when 
including the stellar wind; in Figure [4] we illustrate how our 
addition of the stellar wind affects loop buoyancy. Where 
cool prominences had negative buoyancy out to great heights 
previously (Figure 21 left panel, solid/dotted lines), we now 
see that including the wind means hot loops have negative 
buoyancy at large distances from the stellar surface (Figure 
21 right panel, dashed/dot-dashed lines). 



3.3 Loop Shapes, Maximum Heights 

The critical parameter we are altering in this analysis is 
the external pressure (and thus the pressure balance of the 
loop with its environs). We have improved on previous mod- 
els that assumed a hydrostatic external plasma pressure 
to include the effect of a wind. The height of the source 
surface (the height beyond which all field lines are open) 
is not a well-constrained parameter, as interpretations of 




20 25 5 
Height [R.l 

Figure 4. A comparison of loop buoyancy using the physics of 
JvB05 (left panel) and our present model including a stellar wind 
(right panel). For our fiducial T Tauri parameters, under the pre- 
vious model physics, loops much hotter than their surroundings 
were negatively buoyant only fairly close to the stellar surface, 
within the source surface. 



multi-wavelength observations imply conflicting values. So- 
lar observations show a white light corona which extends to 
greater heights than the X-ray bright corona. In the stellar 
case, densities derived from X-ray emission measures sug- 
gest more compact coronae, while the lack of apparent ro- 
tational modulation in some cases, as well as the formation 
of large-scale prominences, indicates the X-ray coronae are 
extended. Our choice reflects a field more extended than 
main-sequence, solar type stars, but less so than a purely 
dipolar field (based on analysis of the relationship between 
field com plexity and the obs erved coronal X-ray emission 
measure, Ijardine et al.|[2006l ). 

In Figure O we show the shape of hot loops when the 
source surface is placed 3 R* from the stellar surface. Within 
the closed corona, the pressure gradient of the external field 
is directed upward, providing support for the loop against 
the downward force of magnetic tension. This set of solutions 
is effectively identical to those found for hot loops by JvB05; 
we too find that equilibrium solutions can be found for both 
negatively and positively buoyant loops out to the source 
surface. The coolest of the loops plotted (30 MK and 40 
MK) are able to reach heights above the source surface and 
above corotation due to the extra support of the wind. 



4 CONCLUSIONS 

We have addressed the issue of mechanical equilibrium for 
large, hot loops. Since cases have been found in which hot 
loops several stellar radii in height seem to not be disc- 
supported, we propose that these loops reach mechanical 
equilibria beyond the closed coronal magnetic field due to 
the stellar wind. We have shown that, for a generic T Tauri 
star, our equilibrium loop heights are consistent with the in- 
ferred loop lengths for the most powerful flares observed by 
the COUP. We also recover solutions for hot loops within 
the range of heights below the source surface anticipated 
by previous models, suggesting that large and small stable 
post-flare loops can occur under similar conditions. In the 
physics of previous models, which were hydrostatic and did 
not consider the role of stellar winds, loops hotter than the 
quiescent corona were confined to heights within the closed 
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Figure 5. Loop heights as a function of footpoint separation. 
Our loop temperatures probe the parameter space illustrated in 
Figure [l] for loop temperatures from 30-200 MK, we find loops 
with heights ab ove the closed c oronal field (consistent with the 
inferences of iFavata et alj|2005h . The loop temperatures are: 30 
MK (solid line), 40 MK (dotted line), 50 MK (dashed line), 60 
MK (dot-dashed line), and 200 MK (triple dot-dashed line). Grey, 
horizontal lines denote the source surface (extent of closed coro- 
nal field) and corotation radii at 3 R* and 5.2 R*, respectively. 
T e is 2 MK, the plasma /3 is 0.05, and the stellar parameters 
used are as described in Section [3] For all loops shown here, p; is 
~10-20% lower than p e . Footpoint separations beyond 1.05 R*, 
though mathematically possible, are non-physical solutions, as 
we have defined the arcade width to be x = ir/k, where k = 3. 
For a fairly narrow range of loop temperatures, we find a wide 
range of possible maximum loop heights, the coolest loops reach- 
ing just beyond corotation, and the hottest loops confined within 
the closed corona. 



coronal field, needing its support. Our present work indicates 
that, given a stellar wind, hot loops can be found at heights 
well beyond the source surface, and even beyond corotation. 

Young stars are known to be very magnetically active 
in the T Tauri phase, with frequent flaring and high X-ray 
luminosities. The apparent abundance of large scale loops, 
favourable conditions for their formation in TTS coronac 
and winds, and potential for consequence on pre-main se- 
quence stellar (and circumstellar) evolution makes under- 
standing these large-scale loops very important. Having ver- 
ified that such large-scale magnetic structures can indeed 
exist, interesting subsequent work will be to determine how 
the loops' disruption can impact circumstellar disc material 
as well as the star itself. Large-scale disruptions of this na- 
ture could effect stellar rotational evolution and activity as 
well as circumstellar disc evolution and planet formation. 
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